Hydrogen adsorption on boron doped graphene: an ab initio study 
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The electronic and structural properties of (i) boron doped graphene sheets, and (ii) the chemisorp- 
tion processes of hydrogen adatoms on the boron doped graphene sheets have been examined by ab 
initio total energy calculations. In (i) we find that the structural deformations are very localized 
around the boron substitutional sites, and there is an increase of the electronic density of states 
near the Fermi level. Our simulated STM images, for occupied states, indicate the formation of 
a bright (triangular) spots lying on the substitutional boron (center) and nearest neighbor carbon 
(edge) sites. Those STM pictures are attributed to the increase of the density of states within an 
energy interval of 0.5 eV below the Fermi level. For boron concentration of ~2.4%, we find that two 
boron atoms lying on the opposite sites of the same hexagonal ring {B1-B2 configuration) represents 
the energetically most stable configuration, which is in contrast with previous theoretical findings. 
Having determined the energetically most stable configuration for substitutional boron atoms on 
graphene sheets, we next considered the hydrogen adsorption process as a function of the boron 
concentration, (ii). Our calculated binding energies indicate that the C-H bonds are strengthened 
nearby boron substitutional sites. Indeed, the binding energy of hydrogen adatoms forming a dimer 
like structure, aside B1-B2, is higher than the binding energy of an isolated H2 molecule. Since 
the formation of H dimer like structure may represents the initial stage of the hydrogen clustering 
process on graphene sheets, we can infer that the formation of H clusters are quite likely not only on 
the clean graphene sheets, which is in consonance with previous studies [Phys. Lett. 97, 186102], 
but also on the B1-B2 boron doped graphene sheets. However, for low concentration of boron 
atoms, the formation of H dimer structures is not expected to occur nearby a single substitutional 
boron site. In this way, the formation (or not) of H clusters on graphene sheets can be tuned by the 
concentration of substitutional boron atoms. 

PACS numbers; 73.20.Hb;73.20.-r 



I. INTRODUCTION 

Two dimensional crystals of sp^ bonded carbon atoms, 
called graphene, exhibit quite different electronic and 
structural properties compared with their counterpart, 
viz.: graphite (graphene sheets stacked in ABABAB... 
arrangement), and carbon nanotubes (rolled up graphene 
sheets forming nanometer wide cylinders). Indeed, in 
a recent experimental study, Novoselov et al. [1] ob- 
served massless two dimensional (2D) Dirac fermions in 
free standing graphene sheets with high crystal qual- 
ity. In this case, the electron scattering processes are 
suppressed, and thus, the charge carriers can propagate 
freely on the graphene surface. Meanwhile, very recently, 
Meyer at al. observed microscopic corrugations on a sin- 
gle layer suspended graphene sheet, due to the carbon 
atoms displaced out of the graphene plane. Such an ex- 
perimental finding explain the thermodynamic stability 
of those isolated "two dimensional" systems, however, 
"further experimental and theoretical studies are needed 
to clarify the detailed mechanism of the corrugations in 
graphene" [3] . 

On graphene sheets, the sp^ orbitals, forming a bond- 
ing states parallel to the graphene surface, connect the 
carbon atoms arranged in a hexagonal lattice. While 



rise to delocalized tt bonding and tt* anti-bonding states. 
The most of the electronic properties and the chemical 
reactivity of graphene sheets are ruled by those tt and 
TT* orbitals. For instance, the already mentioned mass- 
less 2D Dirac fermions, and the recently observed control 
of the electronic properties in thin films composed by bi- 
layer graphene sheets Q . Focusing on the chemical prop- 
erties, graphene has been considered as a candidate for 
hydrogen storage media [3, Q . Scanning tunneling micro- 
scope (STM) pictures, of hydrogen adatoms on graphene 
surfaces, indicate a local electronic charge enhancement 
and long range electronic perturbation on the graphene 
sheet, giving rise to a local (a/3 x V3)i?30° surface peri- 
odicity Hydrogen on graphene surface has been the 
subject of intense studies, not only addressing the real- 
ization of fuel cells based upon hydrogen storage, but also 
the formation of H2 molecules in the interstellar medium 
(which is an important problem in astrophysics). 

The energetic stability and the equilibrium geome- 
try of hydrogen adatoms on graphene surfaces were in- 
vestigated in detail by Sha and Jackson Q and Du- 
plock et al. @. They find binding energies of 0.67 
and 0.76 eV/H atom, respectively, for a single hydro- 
gen adatom chemisorbed on the graphene sheet. Hav- 
ing characterized the energetic and structural proper- 
ties of hydrogen adatoms on graphene, further experi- 
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molecules through H-H recombination processes on the 
graphene surface. For instance, the Eley-Rideal reaction 
involving free hydrogen atoms (gas phase) and hydro- 
gen adatoms chemisorbed on the graphene surface 0, [9[ . 
On the other hand, the H-H recombination through hy- 
drogen diffusion, on the graphene surface, represents an- 
other reaction path to obtain H2 molecules from the 
chemisorbed hydrogen adatoms ■ IT^. Based upon 

those H-H reaction processes, recent experimental stud- 
ies, supported by ab initio calculations, examined the 
formation of hydrogen clusters on graphene surfaces [l3| . 
Their findings suggest the formation of H dimer like 
structures as an initial stage of the observed hydrogen 
clusters chemisorbed on the graphene surface. The cal- 
culated binding energy of hydrogen adatoms forming a 
H dimer structure is higher compared with the binding 
energy of two isolated H monomers. 

The binding energy of hydrogen adatoms chemisorbed 
on graphene sheets, as well as the formation of H clusters, 
can be tuned by suitable doping processes. For instance, 
the presence of substitutional boron atoms on graphene 
improve the hydrogen adsorption on the carbon atoms 
neighboring the boron substitutional sites [3, [lB| ■ How- 
ever, those calculated hydrogen binding energy results, 
as a function of substitutional boron concentration, are 
quite contradictory. Based upon Raman spectroscopy 
measurements, Endo et al. JT(?| characterized the equilib- 
rium geometry of substitutional boron atoms on graphite. 
For boron concentration of ^^2.7% they find homogeneous 
distribution of substitutional boron atoms. Meanwhile, 
for higher concentration of substitutional boron atoms, 
about 17%, a detailed experimental study, performed by 
Hach et al. [13], indicates the formation of CqB config- 
uration. In the CqB structure, the substitutional boron 
atoms occupy the opposite sites of a carbon hexagonal 
ring. Further experimental/theoretical studies have been 
done aiming to improve our knowledge related to boron 
doped graphite systems [H, [l^, HO] ■ 

In the present paper we report an ab initio total en- 
ergy investigation of boron doped graphene sheets, and 
the chemisorption processes of hydrogen adatoms as a 
function of boron concentration. The presence of sub- 
stitutional boron atoms increases the electronic density 
of states near the Fermi level. Our calculated STM im- 
ages indicate the formation of bright spots on the boron 
substitutional sites, which is in agreement with the ex- 
perimentally obtained STM pictures [3] ■ For higher con- 
centration of substitutional boron atoms, around 2.4%, 
we find that the two boron atoms lying on the opposite 
sites of the same hexagonal ring represents the energeti- 
cally most stable structure. This geometry represents the 
building block of the CqB structure The formation 
of H dimer like structure on graphene has been confirmed, 
as well as the increase of hydrogen binding energy nearby 
the boron substitutional sites. However, in contrast with 
the previous studies, here we have considered the most 
likely configuration for the substitutional boron atoms 



H dimer structure, neighboring a single boron substitu- 
tional site, is energetically less favorable compared with 
two isolated H monomers. Meanwhile, the formation of 
H dimer structure is expected for higher concentration 
of boron atoms, i.e., two substitutional boron atoms oc- 
cupying the opposite sites of a hexagonal carbon ring, 
B1-B2. In this case, the binding energy of H adatoms is 
higher compared with the binding energy of an isolated 
H2 molecule. 



II. METHOD OF CALCULATIONS 

Our calculations were performed in the framework of 
the spin polarized density functional theory (DFT) [2l| . 
within the generalized gradient approximation due to 
Perdew, Burke, and Ernzerhof (2^ . The electron-ion in- 
teraction was treated by using norm-conserving, ab ini- 
tio, fully separable pseudopotentials [2^. The Kohn- 
Sham wave functions were expanded in a combination of 
pseudoatomic numerical orbitals '2^. Double zeta basis 
set including polarization functions (DZP) was employed 
to describe the valence electrons [2^ . The self-consistent 
total charge density was obtained by using the SIESTA 
code Graphene sheets were described within the 

supercell approach, by a single layer of graphene with 
84 atoms, separated by 15 A from their image. A mesh 
cutoff of 170 Ry was used for the reciprocal-space expan- 
sion of the total charge density, and the Brillouin zone 
was sampled by using up to 15 special k points. We have 
verified the convergence of our results with respect to the 
number and choice of the special k points. All atoms of 
graphene sheets were fully relaxed within a force conver- 
gence criterion of 20 meV/A. 

III. RESULTS AND COMMENTS 
A. Boron 

Initially we examined the equilibrium geometry and 
the electronic properties of boron doped graphene sheet. 
For a single substitutional boron atom, Bl in Fig. [Ija), 
we obtained a B-C bond length (dB-c) of 1.50 A, while 
the nearest neighbor C-C bonds are slightly compressed 
(by ^0.01 A) compared with the undoped system. The 
structural deformations of graphene are very localized 
around the substitutional boron site. Experimental STM 
measurements indicate de-c of 1.59 A [l9| which is in 
agreement with our calculated results. 

Graphene sheets are metallic, as indicated by density 
of states (DOS) diagrams shown in Figs. [2^ a) and[2td), 
dashed lines. Such a metallic character has been kept for 
boron doped graphene (solid lines), however, the elec- 
tronic density of states near the valence band maximum 
increases compared with the pristine system (shaded re- 
gions) . Figure (H^b) presents the localization of the elec- 
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FIG. 1: Structural models of boron doped graphene sheet for, 
(a) a single substitutional boron atom (Bl), and (b) two sub- 
stitutional boron atoms per unit cell (51 — 52), corresponding 
to boron concentrations of ~1.2 and 2.4 %, respectively. 



the Fermi level {Ep — 0.5 eV). In this diagram, the 
higher lying occupied states are projected onto a par- 
allel plane at 1.25 A from the graphene sheet. Within 
the TersofF-Hamman approach [27], this projected local 
density of states [Fig. [Hb)[ corresponds to the STM im- 
ages of boron-doped graphene sheet. The bright protru- 
sion on the substitutional boron (center) and the nearest 
neighbor carbon atoms (edges) indicate that the occupied 
electronic states, within Ep — 0.5 eV, are concentrated 
around the boron substitutional site. Those bright spots 
come from tt bonding states of 2^^ orbitals along the B- 
C bonds, see Fig. [21(c) . The formation of bright regions 
around the substitutional B atoms is in accordance with 
the experimentally obtained STM pictures [T9| . 

Within our supercell approach, a single substitutional 
boron atom corresponds to a boron concentration of 
~1.2 %. Two boron atoms per supercell corresponds to a 
boron concentration of ~2.4 %. In this case, we have con- 
sidered several plausible Bi-Bj configurations, indicated 
in Fig. [T](b), for boron atoms occupying substitutional 
sites on graphene sheet. Our total energy results, sum- 
marized in Table HJ reveal the formation of preferential 
domains, or preferential configurations for substitutional 
boron atoms on graphene, viz.: two boron atoms occu- 
pying the opposite sites of the same hexagonal ring, Bl- 
B2 in Fig. [IJb) , represents the energetically most stable 
configuration. Whereas B1-B8 is the energetically least 
favorable structure. The former geometry, B1-B2, was 
proposed by Hach et al. for "boron rich carbon struc- 
tures" [l3|- The structural models Bl-Bi and B1-B5 
are slightly less favorable compared with B1-B2. Since 
in Bl-BA and B1-B5 the boron atoms are placed far 
from each other, d(B-B) — 7.5 and 8.6 A, respectively. 
We can infer that at low concentration regime, the sub- 
stitutional boron atoms may spread out the graphene 
sheet. Whereas, under B-rich condition the formation of 
B1-B2 becomes dominant, giving rise to CqB like struc- 




FIG. 2: Calculated density of states (a)/(d), STM images for 
occupied states within an energy interval of 0.5 eV below the 
Fermi level (b)/(e), and electronic distribution of the highest 
occupied states {Ep — 0.5 eV) along the C-B bonds (c)/(f) 
(unit: e/bohr^), for the structural models BI/BI-B2. 



TABLE I: Total energy differences, AE (eV/B atom), and 
B-B distances in A. 



Configuration 



AE 



d(B-B) 



B1-B2 
B1-B3 
Bl-BA 
B1-B5 
B1-B6 
B1-B7 
B1-B8 



0.00 
0.17 
0.02 
0.01 
0.08 
0.12 
0.66 



2.87 
2.53 
7.49 
8.64 
4.33 
3.86 
1.58 



of the proposed CqB structure. 

Figure[2Id) (solid line) presents our calculated DOS for 
the structural model B1-B2. We observe that the elec- 
tronic density of states increases near valence band max- 
imum (shaded region) compared with the undoped sys- 
tem (dashed line). The occupied electronic states within 
Ep — 0.5 eV are mainly localized on the substitutional 
boron and neighboring the carbon atoms, giving rise to 
the two bright protrusions shown in Fig.[2Kc). Figurc[2Kf) 
indicates that such a STM picture originates from the tt 
orbitals along the B-C bonds. 

It is worth to pointing out that the unlikely B1-B3 and 
B1-B7 fl^ and B1-B8 [3 arrangements {AE = 0.17, 
0.12, and 0.66 eV/B atom, respectively) were considered 
as substrates for hydrogen adsorption on boron doped 
graphene. Thus, we believe that further studies are nec- 
essary addressing the role played by substitutional boron 
impurities for hydrogen adsorption processes, now cou- 
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the substitutional boron atoms, B1-B2. 



B. Hydrogen 

Having established the energetically most probable 
structures for substitutional boron atoms on graphene, 
viz.: structural models Bl [Fig. [Ija)] and B1-B2 
[Fig. mb)], in this section we examined the hydro- 
gen adsorption processes on the clean and boron doped 
graphene sheets. 

Firstly, using the 84-atoms unit cell [Fig.[TJa)], we cal- 
culate the binding energy of hydrogen adatoms on clean 
graphene. For a single hydrogen adatom we find a bind- 
ing energy {E'}j) of 0.98 eV/H-atom. The C-H equihb- 
rium bond length is 1.14 A, which indicates the forma- 
tion of a C-H covalent bond, while the carbon atom un- 
derneath the hydrogen adatom moves upward by 0.5 A 
from the flat graphene layer. For an isolated hydrogen 
molecule, we calculate a binding energy of 4.34 eV/H2- 
molecule (2.17 eV/H-atom). Thus, with respect to the H2 
molecule, the hydrogen adsorption on graphene is an en- 
dothermic process by 1.19 eV/H-atom. Our total energy 
and equilibrium geometry results are in agreement with 
previous ab initio studies of hydrogen adsorption on the 
graphene sheet, for instance, of 0.7-0.8 eV/H-atom 
and C-H bond length of 1.13 A 0,0, [13. 

In a very recent study, Hornekaer et al. observed that 
the presence of a hydrogen adatom on graphene surface 
gives rise to preferential adsorption sites (nearby the al- 
ready adsorbed hydrogen adatom) for the subsequent 
hydrogen sticking process [l^. Indeed we find that a 
hydrogen adatom lying on ha increases the binding en- 
ergy of the next hydrogen adatom (E^^j) on hd or he 
[cf. Fig. [Ha)]. We calculate E'^^h = 2.03 eV/H-atom 
for both hydrogen configurations, ha/hd and ha /he, 
while Horneka;r ct al. obtained binding energies of 1.9 
and 2.1 eV/H-atom, respectively. Those results indicate 
that the formation of ha/hd or ha /he hydrogen dimer 
like structures are energetically more favorable compared 
with two isolated ha monomers, that is, ha+ha — > ha/hd, 
ha /he are exothermic processes. Our findings are in ac- 
cordance with previous studies related to the formation of 
hydrogen clusters on graphene surface. In addition, our 
calculated binding energies E^^jj (for ha/hd and ha /he) 
are comparable with the binding energy of an isolated H2 
molecule (2.17 eV/H-atom). Suggesting that presence of 
a single hydrogen adatom (ha) on the graphene sheet, 
somewhat promote the H2 dissociation nearby ha. 

The binding energy of a single hydrogen adatom, 
-E^, increases upon the presence of substitutional boron 
atoms on the graphene surface. Our results of binding en- 
ergies are summarized in Table HIl For the ha/ Bl config- 
uration, depicted in Fig.[3{a), we find i?^ = 1.89 eV/H- 
atom, whereas on the clean graphene surface we have 
E^fj = 0.98 eV/H-atom. So that, we can infer that 
presence of a substitutional boron atom, in the vicin- 



TABLE IL Binding energies {E'^) of hydrogen adatoms on 
clean and boron doped graphene surface in eV/H-atom. 



Model 
















ha hb 


ha/hb 


ha/hb' 


ha/ he 


ha/hd 


ha/ he 


clean 


0.98 - 








2.03 


2.03 


Bl 


1.89 1.51 


1.36 


1.51 


1.00 


1.67 


1.25 


51-52 


2.13 - 




1.73 




2.37 


1.64 




FIG. 3: Equilibrium geometry and total charge density of 
hydrogen adatoms on the Bl graphene sheet, (a)-(b) ha 
and (c)-(d) ha/hd models. Partial electronic charge density 
within Ef — 0.5 eV for (e)-(f) ha and (g)-(h) ha/hd models. 
Unit, e/bohr^. 



H bond by 0.91 eV. Meanwhile, the hydrogen adsorp- 
tion on top of the substitutional boron atom, hb/Bl, 
is energetically less favorable by 0.38 eV/H-atom com- 
pared with ha/Bl. At the equilibrium geometry, the 
C-H bond length is 1.15 A and the C atom (underneath 
the hydrogen adatom) moves upward by 0.4 A. The for- 
mation of C-H covalent bond is depicted in Fig. El^b). 
Figures ^e) and [3{f) present the occupied electronic 
states along the C-H and B-C bonds, respectively, within 
Ep — 0.5 eV . Those electronic states are composed by 
TT orbitals of carbon and (substitutional) boron atoms 
on the graphene sheet, and spz cr hybridization along the 
C-H bond [Fig.[3Ke)]. Comparing the electronic densities 
depicted in Figs.[2I[c) and[3]Jf), Bl and ha/Bl structures, 
respectively, we can infer that the hydrogen adsorption 
on ha reduces the electronic density of occupied tt or- 
bitals along the B-C bonds. 

Keeping the Bl configuration adsorbed by a hydro- 
gen adatom (ha/Bl), the subsequent hydrogen adsorp- 
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plausible configuration for two hydrogen adatoms on Bl 
graphene sheet. Our binding energy results for the hydro- 
gen adatom on the ha/ Bl graphene sheet, E'^fj, indicate 
that ha /he represents the energetically least stable con- 
figuration. We find E'Ijjj = 1.00 eV/H-atom, nearly the 
same hydrogen binding energy on the clean graphene, i?^ 
= 0.98 eV/H-atom. This result is somewhat expected, 
since he is far from the substitutional boron site, Bl. 
Figure ^c) presents the energetically most stable con- 
figuration, ha/hd, E^fj^ = 1.67 eV/H-atom. The C-H 
bond length is equal to 1.14 A, and its covalent character 
is depicted in Fig. E^d). Similar to the ha/Bl system, 
the electronic states within Ep — 0.5 eV are composed 
by spz <J hybridization along the C-H bonds, and tt or- 
bitals of carbon and boron atoms. Figs. Wig) and[3jh), 
respectively. However, in this case, the adsorption of the 
second hydrogen adatom [hd) is not favored by the pres- 
ence of the first one (/la), E\ju is lower than E\j. Con- 
sequently, ha+ha — > ha/hd is an endothermic process, 
thus, indicating that the formation of hydrogen clusters is 
not expected to occur nearby single boron substitutional 
sites. Thus, for higher coverage of hydrogen adatoms, 
we can infer that hydrogen clusters on the clean region 
graphene surface, i.e. far from boron substitutional sites 
(with E'^jj « 2.0 eV/H-atom) are more likely than hy- 
drogen adatoms neighboring an isolated boron substitu- 
tional site (i^^ff ~ 1.7 eV/H-atom). 




FIG. 4: Equilibrium geometry and total charge density of 
liydrogen adatoms on the B1-B2 graphene sheet, (a)-(b) ha 
and (c)-(d) ha/hd models. Partial electronic charge density 
within Ef — 0.5 eV for (e)-(f) ha and (g)-(h) ha/hb models. 
Unit, e/bohr^. 



Increasing the concentration of substitutional boron 
atoms, forming the energetically most probable config- 
uration, i.e. B1-B2 [Fig. mb)], the binding energy 
of hydrogen adatoms increases compared with the Bl 
graphene surface. For a single hydrogen adatom, ha/Bl- 
B2 in Fig. |Ha), we find equal to 2.13 eV/H-atom, 
which is comparable with the calculated H2 binding en- 
ergy (2.17 eV/H-atom). At the equilibrium geometry, the 
C-H bond length is 1.15 A, and the formation of C-H co- 
valent bond is depicted in Fig. lU^b). Figure |4Ke) shows 
that the spz a orbitals does not lie within Ep — 0.5 eV. 
Whereas the Pz orbital of the carbon atom C(3) strongly 
contributes to the electronic states within this energy in- 
terval. The former cr bonding states, composed by hy- 
drogen s and carbon orbitals, are resonant within the 
valence band of the graphene sheet. In addition, com- 
paring[2If) andjU^f) we verify that electronic density of tt 
orbitals along the B-C bonds increases (asymmetrically) 
upon hydrogen adsorption. 



Different from ha/hd on the Bl graphene sheet, the 
formation of ha/hd dimer like structure on B1-B2, 
namely ha+ha — )■ ha/hd, is an exothermic process by 
0.24 eV/H-atom. Furthermore, E'^^^ is higher than the 
binding energy of H2 molecule, that is, the formation of 
ha/hd on the boron doped B1-B2 graphene surface is an 
exothermic process by 0.20 eV/H-atom compared with 
free H2 molecules. Thus, we can infer that the formation 
of hydrogen clusters can be improved by the increase of 
substitutional boron atoms forming the CqB like struc- 
ture on the graphene surface. Figure [D(d) presents the 
total charge density along the C-H covalent bonds. The 
carbon atoms bonded to hydrogen adatoms exhibit sp'^ 
like hybridizations, where C(3) and C(4) are displaced 
upward by 0.94 A with respect to the pristine graphene 
sheet. The C(3)-C(4) equilibrium bond length (1.54 A) is 
stretched by 0.1 A compared with the C-C bond length 
of undoped systems. The localization of the occupied 
electronic states, within Ep — 0.5 eV, are depicted in 
Figs. [3Kg) andUKh). In the latter diagram we find that 
the electronic density of the tt orbitals along the B-C 
bonds is reduced compared with the tt orbitals of ha/Bl 
[Fig. [4Kf)]. Figure [4jg) depicts the spz cr orbitals along 
the C-H bonds. In a STM measurement, those spz states 
will give rise to bright spots on the hydrogen adatoms ly- 
ing on the B1-B2 boron doped graphene sheets. Indeed, 
recent STM pictures, supported by ab initio simulations, 
indicate the formation of bright protrusions attributed to 
the hydrogen adatoms (forming a dimer like structure) 
on clean graphene sheets [llj . On the other hand, quite 
different STM pictures are expected for ha/ B1-B2. Fig- 
ure Sfe) indicates that, for STM bias voltage of ^0.5 V 
below the Fermi level, the tunneling current from the hy- 
drogen adatoms will be negligible compared with the one 
from neighbor carbon and boron atoms. Thus, in this 

case, the hydrogen adatom site will appear darker than 
1 1 1 1 _j 
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IV. CONCLUSIONS 

In summary, we have performed an ab initio total en- 
ergy investigation of boron doped graphene sheets, and 
their interaction with hydrogen adatoms. For a single 
substitutional boron atom, the equilibrium geometry of 
graphene sheet is weakly perturbed, while the electronic 
density of states near the valence band maximum has 
been increased. Our simulated STM images indicate the 
formation of bright (triangular) spots on the boron sub- 
stitutional site (center) and the nearest neighbor carbon 
atoms (edges). For two boron atoms occupying the sub- 
stitutional sites on the graphene sheet, corresponding 
to a boron concentration of ~2.4%, we find that boron 
atoms lying on the opposite sites of the same hexago- 
nal ring {BI-B2 structure) represents the energetically 
most stable configuration. Furthermore, the hydrogen 
adsorption on graphene has been examined as a func- 
tion of the concentration of boron atoms. For the pris- 
tine graphene system, wc find that the binding energy of 
hydrogen adatoms forming H dimer like configurations 
(ha/hd or ha/he) is higher by 1 eV/H-atom compared 
with the binding energy of two isolated H monomers {E'^ 
= 0.98 eV/H-atom). This is in accordance with the re- 
cent experimental/theoretical study of the formation of 
H clusters on graphene sheets [ll|, [l^ ■ The presence of a 
single substitutional boron atom increases the hydrogen 
binding energy by 0.9 eV/H-atom compared with the one 
on the clean graphene sheet. In this case, the hydrogen 



adatom lies on top of the carbon atom nearest neighbor 
to the substitutional B site [ha/Bl configuration). How- 
ever, different from the pristine system, the adsorption 
of two hydrogen adatoms close to each other, forming 
a H dimer like structure, is energetically less favorable 
compared with two isolated H monomers. This result 
indicates that the formation of hydrogen clusters on the 
graphene sheet is somewhat suppressed by the presence 
of a single substitutional boron atom. Finally, for boron 
concentration of ^2A% we find that the binding energy 
of H dimer structure {E\jjj) is higher than the binding 
energy of an isolated II2 molecule, where the hydrogen 
adatoms lie on the carbon atoms neighboring the substi- 
tutional boron sites. In this case, boron atoms give rise 
to preferential adsorption sites for hydrogen adatoms on 
graphene sheets, and thus, improving the formation of 
hydrogen clusters on graphene sheets. Those results indi- 
cate that the formation (or not) of H clusters on graphene 
sheets can be tuned by the concentration of substitu- 
tional boron atoms. 
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